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Preface 



This report is a collection of material that has been used 
in courses on search, detection and localization modeling. Its 
organization follows to some extent material by S. M. Pollock in 
Selected Methods and Models in Military Operations Research which 
is listed in the report bibliography. The report is not intended 
to be a text on these subjects. In particular, in some areas it 
does not provide the depth of coverage that is found in the book 
Search and Detection by Alan R. Washburn which is cited in this 
report as Reference 22. 

In the third revision, typographical and other errors have 
been corrected and, in addition, changes and additions have been 
made to Sections XII and XIV of the report. 
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I. Detection Models and Signal Detection Theory 

Signal detection theory is the basis for analyzing the 
detection models that are described in this report. In signal 
detection theory, the decision making portion of a detection 
system is called the receiver and a detection experiment is the 
observation by a receiver of input data accumulated during some 
time interval. The data that is related to a target is called 
signal. The data that is not related to the target is called 
noise. In general, the target data is associated with a 
localization region that in some cases is called a resolution 
cell. When a detection experiment is performed, either the event 
Hq = {the receiver input is noise} or its complement 

= (the receiver input is signal and noise} will occur. In 
the first detection models that are described here, after 
analysis of the input data by a receiver, either the event 
Dq ={the receiver decides the input is noise} or its complement 
D]^ = (the receiver decides the input is signal and noise} will 
also occur. Detection models for which Dj^ is the complement of 
Dq are called binary detection models or forced choice detection 
models. Four events which are important in binary detection 
models are indicated in the Venn diagram of Figure 1. 

The Venn diagram emphasizes a decision problem that is 
associated with a receiver that can be modeled using a binary 
detection model. The problem is this: Under what conditions 

should the event occur? That is, under what conditions 
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should a receiver decide that the input data accumulated during 
the observation time interval is signal and noise? 



Ho 



Hi 



Dq Hq 


Dq n 


Di n Hq 


□i n Hi 



Figure 1. Four events of importance in binary detection models. 

In the detection model descriptions that follow, the 
following notation and terminology are used: pf = P(Di|Hq), the 

probability of given Hq, is called the false alarm 

probability; p^ = P(D2|H2), the probability of given 

is called the detection probability and P = P(H]^), the 
probability of is called the prior probability. 

In the detection models, the input to a receiver is 
determined by a stochastic process that has the following 
characteristics: It is a random noise process when there is no 

target data and it is a random noise process plus a signal 
process when there is target data. Although the receiver input 
process in some cases might appear to be determined by a 
continuous parameter stochastic process, because of the finite 
amount of information (unitjue data) contained in a bounded 
sequence of finite length, a discrete parameter stochastic 
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process is sufficient to determine the receiver input in these 
cases. This is established formally by the stochastic sampling 
theorem. Consequently, in these models, the input to a receiver 
is determined by a sequence of random variables Yj^, ••• ,Yjj and 
an observation yields a sequence of values y^, . . - Yin* 

Three detection models are described in Section III. In the 
first model, the signal process is a deterministic process. In 
the second and third, the signal process is a random process. To 
define a random noise process or a random signal process, only 
the joint distribution of the finite sequence of random variables 
that determine the process needs to be specified. If the signal 
process is a deterministic process, the signal values can be 
determined before an observation is performed. To define the 
process in this case, only these values need to be specified. 
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II. Decision Criteria 

To simplify the discussion of decision criteria and 
decision rules, a receiver's input will be assumed to be 
determined by a single decision random variable Y. In this 
case, the input process in determined by the conditional 
distribution function Fy(y|Ho) when the input is noise alone 
and by the conditional distribution function FyCylHi) when the 
input is signal plus noise. 

The condition that a receiver's input is required to 
satisfy in order that the event will occur can be specified 

in terms of a decision rule. For the assumed case, a decision 
rule is a rule which determines for every observable value of Y 
the decision that the receiver is to make. The decision rule can 
be considered to be a function 0(y) which relates each 
observable value y to one or the other of the following two 
commands: do = "decide that the receiver input was noise" and 

d^ = "decide that the receiver input was signal and noise". 
Choosing a decision rule 0(y) defines a set n of observable 
values of Y such that the event = { Y e n } . 

The problem which was considered in Section I can now be 
restated in the following way: What criterion should be adopted 

in order to determine a decision rule or, equivalently, its 
corresponding set n ? A desirable characteristic for a criterion 
is suggested by the following argument: Consider the odds in 

favor of given y is observed. That is, consider 
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P(Hi|Y = y)/P(Ho|Y = y) . One might expect that y would be a 
member of the set n if and only if y made this ratio equal to 
or greater than some value k. But this is equivalent to defining 
n as follows: n = { y ; L(y) ^ K ) where L(y) is the 
likelihood ratio associated with an observed value y and K is 
a constant related to the constant k. This suggests that 
choosing an optimum criterion is equivalent to choosing an 
optimum value for K. Four specific decision criteria are defined 
next in terms of K. For each criterion, n has the above 
form. But for each criterion the choice of K is different. 

The decision criteria are: 

1. The Neyman-Pearson Criterion: Choose n so that p^j 

is a maximum subject to the constraint that pf < a where a is 
a specified value. For a continuous decision random variable, 
the constant K is chosen so that pf = a. 

2. The Bayes Criterion: Choose fl so that the expected 

cost of a receiver's decision is a minimum. For a continuous 
decision random variable, K = t (Cio~‘^00)/(‘^01~‘^ll) ] * 

Cio > Cqo and cqi > where c^j is the cost of n H j . 

3. The Ideal Observer Criterion: Choose n so that the 

probability that the receiver makes an incorrect decision is a 
minimum. K = (1-P)/P for a continuous decision random variable. 

4. The Minimax Criterion: Choose n when P is unknown so 

that the maximum expected cost of a receiver's decision is a 
minimum. If cj^q > <^00 ^01 ^ *^11' then 
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K = [ (cj^o'^^oo)/ ^ continuous decision 
random variable. Here, P* is the value of the prior probability 
P that would make the expected cost of a receiver's decision a 
maximum if P were known and the Bayes Criterion were used. 

If a model which specifies the conditional distributions 
Fy(y|Ho) and FY(y|Hj^) and a decision jrule are adopted, then 
the value of pf and the value of p^j are determined. This 
pair of values (Pf,p<j) is called a receiver operating point. 

If the decision rule results from using a likelihood ratio 
criterion such as one of the four listed above, then it will 
involve the parameter K since n = (y: L(y) > K}. And, for a 
given value of K, since n uniquely determines the pair 
(Pf,Pd) , a single operating point results. By varying K, a set 
of operating points can be generated which determines a receiver 
operating characteristic curve or ROC curve. Different ROC 
curves can be produced by changing either one or both of the 
conditional distributions which implies either a change in the 
signal process or a change in the noise process. 

A decision rule which results from using a likelihood ratio 
criterion in a model in which the input process is determined by 
a set of m random variables can be expressed in terms of a set 
n as follows; n = { (y^, ••• ,yj,) : L(yi, ••• ,yjn) > K ) where 
K is specified in the same way that it is when m = 1. 
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III. Three Binary Detection Models 

Three detection models are examined in this section. 

For the first two detection models, the input stochastic process 
for an observation is defined by a time sequence of continuous 
random variables. The random variables represent a sample from a 
continuous parameter stochastic process which is sampled at times 
such that the random variables are independent. For the third 
detection model, the input stochastic process is a counting 
process and it is defined by a single discrete random variable 
that is equal to the number of events that are counted during the 
observation. 

Model I: In the first detection model, a sampled noise 

value is a value of a normally distributed random variable with 
mean zero and with known variance a*. And a sampled signal 
value is a known value of a deterministic variable. Thus, the 
input process corresponding to an observation consists of some 
number m of independent normal random variables Y^, ••• ,Yjn 
each with variance o*. And, for i=l,2, ••• ,m, when a signal 
is not present the mean of is zero and when a signal is 

present the mean is sj^. The result of using a likelihood ratio 
decision rule in the model can be expressed in terms of a random 
variable Z. This random variable is called a crosscorrelation 
statistic and it is defined by Z = S s^'Y^ where the sum index 
1 = 1,2, ••• ,m here and in the remainder of this section. 
However, it is more convenient to express the result in terms of 
a statistic V which is defined by V = Z/a 2 . In terms of this 
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random variable, the conditional probabilities pf and P(j are 
given by: pf = 1 - *(v*) and p^j = 1 - *(v* - d^) where ♦ 

symbolizes the standard normal cumulative distribution function, 

V* = (l/a 2 )*(a*ln K + (1/2) E , is determined by the decision 

2 2 , , . 
rule and d = S Sj^/o is called the detection index. 

Often, the input stochastic process represents a quantity 
whose square is proportional to power. In such a case, the 
average receiver input power is the random variable E Y^/m. If 
a signal is not present, the expected average receiver input 

power is N=Eo*/m=o* where N is called the noise. The 

. . . . * 
average receiver input signal power is S = E s^/m where S is 

called the signal. In these terms, d = m* (S/N) where S/N is 

called the signal-to-noise ratio. 

If a receiver's input can be considered to be a time 
sequence of continuous voltage values such as in the case of a 
sonar receiver, in some cases a frequency representation can be 
used that involves the concept of receiver bandwidth. In these 
cases, the noise process is assumed to be such that m = t/5t 
where t is the integration time (the duration of an 
obseirvation) and St is the time between samples with 
St = 1/[2(BW)] where BW is the bandwidth and St is 
determined by the sampling theorem. This implies that the 
detection index can be written as d = 2t* (BW) • (S/N) . By 
defining Nq as the power spectral density where Nq = N/BW^ the 
detection index can also be written as d = 2t* (S/Nq) • 
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In Reference 2, the conditions required by the first model 
are called Case I and in the following sections the first model 
is called the Case I model. A receiver that processes data 
such that it would implement a likelihood ratio decision rule 
under the conditions of the first model is called a matched 
filter or crosscorrelation detector. If the description of the 
input noise is adequate, a Case I model can be used to obtain an 
estimate of an upper bound on a detection system's performance, 
since all the information necessary to define the signal is 
assumed to be known. 

Model II: In the second detection model, a sampled noise 

value is a value of an independent normal random variable with 

mean zero and with known variance a*. And a sampled signal 

value is an independent random variable with mean zero and with 

known variance a|. Thus, the input process corresponding to an 

obseirvation consists of some number m of independent normal 

random variables •••,’/!„ each with mean zero and each with 

variance a* when a signal is not present and each with variance 

a* + Og when a signal is present. The result of applying a 

likelihood ratio decision rule in this model can be expressed in 

2 

terms of a statistic X which is defined by X = S Y^. 

When a signal is not present, the statistic X/N has a 
chi-square distribution with m degrees of freedom. When a 
signal is present, the statistic X/(N+S) has a chi-square 
distribution with m degrees of freedom. So, in terms of these 
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two statistics, the two conditional probabilities pf and P(j 
are: pf = P(X^ > x*/N) and p<j = P{X^ > (x*/N) • [ 1/ ( 1+S/N) ] ) 

where is a chi-square random variable with m degrees of 

freedom, x* is a number which is determined by the decision rule 
and S/N is the signal-to-noise ratio. A receiver that would 
implement a likelihood ratio decision rule under the conditions 
of the second model is called an energy detector or square law 
detector. 

The mean of a chi-square random variable with m degrees of 
freedom is m and the variance is 2m. By the central limit 
theorem, as the number of degrees of freedom of a chi-square 
random variable becomes large, it can be approximated by a normal 
random variable that has the same mean and variance. Therefore, 
Pf and P(j can be approximated by pf = 1 - §[(x*/N - m)/(2m)i] 
and Pd = 1 - *{ [l/( 1+S/N) ] [x*/N - m - m* (S/N) ]/ (2m) H • And, 
with V* = [x*/N - m)/(2m)h and d = (m/2) • (S/N) * , this 
becomes; pf = 1 - *(v*) and p^ = 1 ” [l/( 1+S/N) ]• (v* - d^ ) . 

Based on the above approximation, if the noise N is 
significantly larger than the signal S, then Pf and P(j can 
be approximated by: pf = 1 - §(v*) and p^ = 1 - f(v* - d^). 
And, if the concept of bandwidth is applicable so that the sample 
size m = 2t* (BW) , then d = t • (BW) (S/N) * . In Reference 2, the 
conditions required for this approximation are called Case II and 
in the following sections this limiting form of the second model 
is called the Case II model. 
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Model III: In the third detection model, a sampled noise 

value and a sampled signal value are values of independent random 
variables that are determined by independent Poisson processes 
that are observed for a time interval t. The noise process is 
characterized by a counting rate a, the signal process is 
characterized by a counting rate and the noise and signal 

processes are additive. This implies that when the input is 
noise alone, the input is a Poisson random variable with 
parameter a*t, the expected number of noise counts, and when the 
input is signal and noise, the input is a Poisson random variable 
with parameter (a + Ug) *t, the expected number of noise and 
signal counts. 

For a likelihood ratio decision rule, pf = 1 - P(y*;a*t) 
and P(j = 1 - P[y*;(a + Og) ‘t] where y* is a threshold value 
that is determined by the decision rule and P(y;©) represents 
the Poisson cumulative distribution function with parameter 0. 
When 6 is large, the cumulative distribution function can be 
approximated by the cumulative distribution function of a normal 
random variable that has the same mean and variance. Using this 
approximation for cases where a*t is sufficiently large, since 
both the mean and variance of a Poisson random variable are equal 
to e, Pf = 1 - §(v*) and Pd = 1 - ♦{ [ 1/ (1 + ag/a)^](v* - d^ ) 
where v* = (y* - at)/(at)^ and d = a*t*(ag/a)*. If, in 
addition, a is significantly larger than Og, that is, if 
Qg/a « 1 as well as a*t » 1 , then pf and Pd can be 
approximated by: pf = 1 - #(v*) and Pd = 1 - *(v* - d^). 
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The third detection model might be used to describe a 
receiver whose input for an observation is the number of photons 
counted by a radiation detector in situations where a*t, the 
expected number of counts when no signal is present, is of the 
order of thirty or more. 

When a likelihood ratio decision rule is used in the three 
models discussed above, for the first model and under limiting 
conditions for the second and third models, the following result 
is obtained: pf = 1 - #(v*) and Pd = i ~ #(v* - d^) where the 
value of V* depends on the noise power N for the first and 
second models. For a sonar receiver described by the first 
model, that is, by the Case I model; d = 2t • (BW) (S/N) . For a 
sonar receiver described under the limiting conditions for the 
second model, that is, by the Case II model, d = t • (BW) (S/N) * . 
So, in either a Case I model or a Case II model of a sonar 
receiver, the detection index d is a function of the time 
bandwidth product t*(BW) and the signal-to-noise ratio S/N. 
Since sonar equations relate S/N to system, target and 
environmental parameters, a sonar equation can be used to relate 
S/N to these parameters in a model of a sonar receiver. 
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IV. General Detection Models 



The detection models that have been considered to this 
point are based on binary detection theory. After each 
observation, a receiver decides either that the input 
corresponding to the observation was noise or else it decides it 
was signal plus noise. However, in some detection systems this 
decision is delayed. In a computational sense, a model of such a 
detection system is generally complex relative to a binary 
detection model. To illustrate this, consider an active sonar 
system whose receiver includes an operator. Suppose the 
probability that the operator will detect a target echo has been 
determined in a laboratory experiment in which the operator was 
required to decide after each input corresponding to a resolution 
cell that either the input was a target echo (signal) and noise 
or the input was noise alone. In addition, suppose that under 
operational conditions the operator normally delays this 
decision. Then, in general, the probability that the operator 
will decide that the input corresponding to a resolution cell 
that contains a target is a target echo and noise will not be 
equal to the probability of the event in the forced choice 
experiment. And, in addition, the probability that the operator 
will decide the input corresponding to a resolution cell that 
does not contain a target is a target echo and noise will not be 
equal to the probability of this event in the forced choice 
experiment. Consequently, in general, the value of both p^ and 
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Pf for an operational environment will be different than that 
for the laboratory environment. 

One model that has been proposed to deal with this kind of 
situation defines the event that a receiver decides that the 
input corresponding to a resolution cell is signal and noise to 
be equivalent to the event that out of n consecutive 
observations at least k of them would result in the decision 
that the input was signal and noise in a forced choice 
experiment. The model is said to be based on an k-out-of-n 
detection criterion. With this criterion, the probability that a 
target will be first detected on the observation can be 
found as follows: Determine the 2^ sequences of forced choice 

responses that could result for a sequence of j consecutive 
observations. Next, determine the probability of occurrence for 
each sequence that first satisfies the k-out-of-n detection 
criterion on the observation. The probability of first 

detection on the observation is equal to the sum of these 

probabilities. The cumulative probability of detection at the 
jth observation is the sum of the probabilities of first 
detection on the i^^ observation for i=l,2, ••• ,j. 
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V. Signal-to-Noise Ratio Detection Models 

In some radar and sonar detection models, for a 
specified value of Pf, a minimum acceptable value of P(j is 
defined. This minimum acceptable value of p^j and the specified 
value of Pf define what can be called a minimum acceptable 
signal-to-noise ratio (S/N)j[, if P(j is a nondecreasing 
function of signal-to-noise ratio. In some sonar detection 
models, (S/N)j,j in decibels is called the detection threshold 
DT. In symbols, DT = 10 log(S/N)m. If the minimum acceptable 
value of P(j is .5, then DT is usually called the recognition 
differential RD. The difference between the signal-to-noise 
ratio in decibels and RD (or DT) is called the signal excess 
SE. In symbols, SE = 10 log(S/N) - RD. 

One interpretation of signal excess is that for a 
localization region containing a target detection occurs with 
probability one if SE > 0 and with probability zero if SE < 0. 
This interpretation provides the basis for defining detection in 
the three encounter detection models that are discussed in 
Section VII. A more consistent interpretation of signal excess 
is: If SE > 0, then the probability of detection p^j is 

greater than or equal to the minimum acceptable value (.5 if 
recognition differential RD is used to define signal excess) . 
For cases where p^j increases rapidly with signal excess in the 
neighborhood of zero signal excess, the two interpretations may 
be operationally equivalent. For a discussion of this point as 
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well as a discussion of an operational case in which receiver 
decisions are delayed, see Reference 3. 

Signal excess (signal-to-noise ratio) detection models 
provide a basis for general detection models, in particular, 
models that describe nonstationary noise and signal processes and 
randomly changing decision rules. This is illustrated by the 
models described in Section VII. In addition, signal excess 
models provide a basis for delayed receiver decision models. 

This is illustrated by the active sonar detection models in both 
Reference 4 and Reference 5 that are based on a k-out-of-n 
detection criterion. In all of these models, the signal-to-noise 
ratio and the recognition differential are random variables. 

Using X(t) to represent a random variable corresponding to 
an index time t and a subscript to identify the random variable 
in such models, for a passive sonar receiver, the signal-to-noise 
ratio in decibels associated with a decision at the index time 
is: XgL(t) - X«pL(t) - [XjjL(t) - Xpj(t)]. In this expression, 

SL represents source level, TL represents transmission loss, 

NL represents noise level and DI represents directivity index. 
Since signal excess SE is defined to be the difference in 
decibels between the signal-to-noise ratio and the recognition 
differential (or detection threshold) , it too is a random 
variable and, for any decision time t, one can write: 

(1) XsE(t) = XsL(t) - XTL(t) - [XNL(t) - XojCt)] - XRD(t). 

The distributions of the random variables on the right side of 
Equation 1 determine the distribution of the signal excess. In 
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the passive sonar detection model described in Reference 6, 

' Xp[)(t) and, in effect, %l(^) normally 

distributed random variables while XrpL(t) is a uniformly 
distributed random variable. In the three signal excess models 
that are described in Section VII, all of the random variables in 
Equation 1 are normally distributed. 

It is sometimes convenient to write Equation 1 as follows: 

(2) XsE(t) = SE(t) + X(t) . 

In Equation 2, SE(t) is the expected value of the signal excess 
determined by the following expected value equation; 

(3) SE(t) = SL(t) - TL(t) - [NL(t) - DI(t)] - RD(t) 

where each term on the right represents the expected value of the 

indicated random variable and X(t) is a random variable that 

determines the stochastic character of the signal excess. Since 

SE(t) is the mean of XgE(t); by Equation 2, the mean of X(t) 

is equal to zero and the standard deviation of X(t) is equal to 

the standard deviation of Xg£(t) . If a represents the 

standard deviation of XsE(t) the random variables on the 

right side of Equation 1 are statistically independent, then 
2 2 2 2 2 2 

a = OgL + cjipL + *^01 ‘^RD* This relation has been used to 

determine a standard deviation for the signal excess in 
operational models. 
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VI. General Encounter Models 



A basic problem associated with search modeling is that 
of determining the probability that a target will be detected by 
a detection system during an encounter with one or more detection 
systems. In the encounter models that are considered in this 
report, during a search, observations are made of a series of 
localization regions. The probability of detection on an 
observation is P(Dj^ n . The probability of a false alarm on 
an observation is n Hq) . In these models, the time to 

resolve a false alai*m is ignored. However, p^ and Pf are 
assumed to be determined by some criterion such that Pf is an 
operationally reasonable value. 

Using the order number of a decision rather than its time as 
an index relative to detection decisions for localization regions 
that contains a target and a random variable N to represent the 
decision order number at which detection first occurs, the 
probability of detection during an encounter can be written as: 
P(N < n) = P(N < m) + P(N = m+1) + ••• + P(N = n) or equally as 

P(N < n) = 1 - [1 - P(N < m)]*(l - g^+i) ••• (1 - gn) where 

= P(N = i|N < i-1) is the probability of the event detection 

at the i^^ decision conditioned on the event no detection at an 

earlier decision and 1 < m < n. The second expression is 
generally of greater interest than the first expression, since 
g^ can usually be more directly related to operational 
parameters such as range and environmental conditions that 
determine a target's detectability than can P(N = i) . 
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with a time rather than the order number to index a decision 
and a random variable T to represent the time index at which 
detection first occurs, P(N < n) becomes P(T < tj^) with 
P(T < t„) = 1 - [1 - P(T < t„,)].[l - g(t^+i)] ••• [1 - g(tn)] 
where g(ti) = P(T = t^jT <ti_i) . 

If g(tj^) « 1 for i = 1,2, ••• ,n, then, to a first 
approximation, ln[l - g(t£)] = -g(t£) for 1 = 1,2, ••• ,n and 
P(T < tj^) = 1 - [1 - P(T < tj„)]»exp[-Z g(ti)]. This follows 
since P(T < t^) = 1 - [1 - P(T < t^) ] •exp[Z ln[l - g(ti) ] where 
the sum index i = m+i, ••• ,n. A continuous analog to this 
approximation can be used to describe an encounter for which 
g(ti) « 1 for i = m,m+l, ••• ,n and decisions during the 
encounter can be considered to occur continuously. That is, the 
time of an observation corresponding to a decision and the time 
between decisions are both negligible relative to the time of the 
encounter. 

The analog can be developed as follows: First, let St be 

the time between decisions, then t^ = i*6t and the probability 
of detection P(T < t^,) = 1 - [1 - P(T < tn,)]»exp[-Z T(ti)*St] 
where T(tj^) = (l/6t) *g(tj^) is a detection rate function 
(a probability of detection per unit time) and, in terms of St, 
the probability g(tj^) = P[T = i*St|T < (i-l)»St]. 

If T is considered to be a continuous random variable, the 
expression for P(T < t^) above indicates that the sum in the 
exponent should be replaced by an integral whose integrand is a 
continuous function T(t). If T(t) can be determined, then. 
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with g(ti) as a guide, the cumulative probability of detection 
P(T < t) can be defined by: 



(4) 



T(t) = lira {(1/St) P(t < T < t+5t|T < t) ) 



where the limit is for St approaching zero. Equation 4 
implies the differential equation: dp(t)/dt = [1 - p(t)]*T(t) 

where p(t) = P(T < t) . A solution to this equation is; 



where t is the time index for a decision during an encounter, tm 
is some time during the encounter and tn > tm. A T(t) that is 
based on a visual detection model is described in Reference 7. 

If the detection capability of a detection system is assumed to 
depend on a target's position relative to the detection system 
during an encounter but not to depend on the clock time, then the 
time index of a decision can be a relative index that determines 
the target position that is associated with a decision rather 
than the clock time associated with the decision. 

The above results apply to the case of an encounter between 
a target and a collection of detection systems. However, if the 
detection systems are not collocated, it is generally convenient 
to describe encounters of this kind in terms of encounters 
between the target and the individual detection systems. In 
either case, if the event target detection for a detection system 
is not independent of the event for other detection systems, then 
in order to describe this in an encounter model the correlation 
between the input to the detection system and the inputs to the 
other detection systems must be specified. This has been done in 



(5) 



P(T < tn) = 1 - [1 - P(T < tm) ] 
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some models as follows: First determine the probability of 

detection for each system acting alone. Let Pj^ be the 
probability that the i^^ system detects the target during the 
encounter under this condition. Next, consider two cases; In 
the first case, the random factors that determine detection for a 
system are independent of those that determine detection for the 
remaining systems. In the second case, the random factors that 
determine detection for the systems are completely dependent. In 
the first case, the probability that at least one system detects 
the target is given by: Pj = 1 - (1 - • (1 - P 2 ) (1 “ Pn) 

where n is the number of detection systems involved. In the 
second case, the probability that none of the systems detect the 
target is given by: 1 - Pp = 1 - Pj„ where Pj„ > Pj^ for 
i = 1,2, ••• ,n since if the m^^ system does not detect the 
target, none of the remaining systems will detect it. The 
probability that at least one system detects the target is given 
by: P= a*Pp+ (1 - a)*Pj where a determines the degree of 
correlation and 0 < a < 1. A way to determine a value for a 
is described in Reference 8 . 
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VII. Three Signal Excess Encounter Models 

In the three models described in this section, 
detection is defined in terms of signal excess as it is in 
Section V. Each model determines a cumulative probability of 
detection for a target in an encounter with a passive sonar 
system. An observation in the models is indexed by time and the 
index can usually be considered to be the time at the end of the 
observation. During an encounter, observations are made of one 
or a series of localization regions. By implication, a false 
alarm can occur for a localization region that does not contain a 
target during an observation since the value of RD (or DT) is 
determined by some specified false alarm probability. However, 
as they are generally used, signal excess models do not account 
for false alarms. This can be viewed as equivalent to modeling 
the time to resolve a false alarm to be effectively zero. 

To determine signal excess in the models, it is convenient 
to use Equation 2. For each decision in an encounter, there is a 
random variable X(t) defined by Equation 2 that determines the 
random character of the signal excess. For a sequence of 
decisions, the set of these random variables ordered by their 
time index constitutes a stochastic process. And the joint 
distributions of these random variables determines the nature of 
the stochastic process. In the three encounter models described 
in this section, the stochastic process is called a lambda-sigma 
jump process. The time series that are generated by lambda-sigma 
jump processes are represented by the plot in Figure 2 below. 
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The jumps in the time series occur at times determined by a 
Poisson process with a mean jump rate lambda. This implies that 
the time between jumps is a random variable with an exponential 
distribution and that the expected times between jumps t is 
equal to the reciprocal of lambda. 




Figure 2. A time series representing a realization of a lambda- 
sigma jump process. On the plot, a in dB equals one unit on 
the vertical axis and t equals one time unit on the horizontal 
axis. 



From Figure 2 , note that the observed values of neighboring 
random variables are equal unless a jump has occurred between 
the observations. When a jump occurs, the first random variable 
after the jump is normally distributed with mean zero and 
variance a* and it is independent of all the random variables 
before the jump. Conditioned on a jump pattern, this random 
variable and all the random variables between it and the next 
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jump are dependent and the correlation coefficient between any 
pair is one. That is, if the value of the signal excess is known 
at some time, then all of the values between the last jump before 
that time and the first jump after that time are also known. 
However, since the jumps occur randomly, knowing the value of the 
signal excess with certainty at some time does not determine the 
values of the signal excess with certainty at neighboring times. 
In the unconditioned case, the correlation coefficient between 
the random variables X(t) and X(t+T) is equal to 1/e. For 
this reason, t is referred to as a relaxation time. 

It appears that the use of the lambda-sigma jump process is 
based more on past practice than on experimental justification. 

In this regard, see Reference 9. By referring to Equation 1, it 
can be seen that the lambda-sigma jump process is determined by 
the sum of the stochastic processes that determine the random 
variables on the right side of this equation. Although the sum 
of a collection of normal random variables is a normal random 
variable, the sum of a collection of lambda-sigma jump processes 
is not a lambda-sigma jump process. This suggests that if the 
lambda-sigma jump process does adequately describe the 
variability of the signal excess, then the majority of the 
variability of the signal excess may be due to a single one of 
its elements. For example, transmission loss. 

In the three encounter models described below, detection is 
defined in terms of signal excess as described in Section IV and 
decisions are indexed by a time that can usually be considered to 
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be the time of the decision. During an encounter, observations 
are made sequentially of one or a series of localization regions 
(resolution cells) . For a localization region that does not 
contain a target, the signal observed during the observation of 
the region is zero. For these observations, the time to resolve 
a false alarm is zero. However, since the value of RD (or DT) 
is finite and consequently the false alarm probability is not 
zero, by implication, the cost associated with a false alarm is 
not zero. 

The First Passive Sonar Encounter Detection Model: This 

model describes an encounter in terms of a series of decisions 
with each decision based on the signal excess X 3 g(t) at a time 
corresponding to the end of an observation. The observations are 
of equal duration and the integration time that determines the 
recognition differential is equal to the duration of the 
observations. In the model, XsE(t) is determined by a lambda- 
sigma jump process. For an encounter of m observations in 
which SE(t) is unimodal and in which the time of the single 
maximum is prior to or at the end of the encounter, it is shown 
in Reference 10 that the probability p that detection will 
occur during the encounter is given by the following equation: 

(6) p = 1 - [(1 - Pc)/(1 - P-Pc)]-(1 - P-Pi) ••• (1 - P*Pm) 
where ^ = 1 - exp(-St/T) and p^ = $[SE(ti)/a] for i = 1,2, 
••• ,m. Here, 6t indicates the duration of an observation and 
* indicates the standard normal cumulative distribution function 
as before. The integer c is the index of a decision time t^ 
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for which SE(tc) is greater than or equal to SE(tj^) for any 
time t^ and t^ < t^ < tj„. 

As T approaches zero, p approaches one and Equation 6 
approaches this form; 

(7) p = 1 - (1 - Pi) • • • (1 - pj„) . 

In this limit, the signal excess random variables are all 
independent. Note that Equation 7 applies without the condition 
that SE(t) be unimodal. 

As T approaches infinity, p approaches zero and 
Equation 6 approaches this form: 

(8) p = Pc- 

In this limit, the correlation coefficient between any pair of 
signal excess random variables is equal to one. Note that 
Equation 8 applies without the condition that SE(t) be 
unimodal. Equation 8 defines a complete dependence encounter 
model . 

The Second Passive Sonar Encounter Detection Model: This 

model is in a sense a third limiting form of the first passive 
sonar encounter detection model. In this limit, the time between 
decisions approaches zero. However, in this limit the 
integration time that determines the recognition differential is 
not equal to 5t and it does not approach zero. It is, in 
effect, chosen by the user of the model through the user's choice 
of the value for the recognition differential. For an encounter 
that begins at tl and ends at t2 and for which XgE('t) is 
determined by a lambda-sigma jump process and SE(t) is 
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unimodal, it is shown in Reference 10 that for this limit, 
Equation 6 has the following form: 



(9) p = 1 - [1 - p(tc) ] *exp[-(l/T) • _ p(t) dt] 

where p(t) = #[SE(t)/a] and where now t^, is the encounter 
time such that SE(t£,) is greater than equal to SE(t) for any 
other encounter time t and tl < t^ < t2 . 

The Third Passive Sonar Encounter Detection Model: This 

model describes an encounter between a target and a passive sonar 
detection system in which detection occurs during an encounter if 
the average value of the square of the continuously observed 
signal-to-noise ratio over a time interval of length u is 
greater than or equal to the square of the signal-to-noise ratio 



that determines the recognition differential for an integration 
time equal to u. With R(s) the random signal-to-noise ratio 
at a time s and the random signal-to-noise ratio that 

determines the random recognition differential for an integration 
time u, detection during an encounter occurs at the first time 
t that the following inequality is satisfied: 



where the time origin is chosen so that t > 0 and where the 
integration time u = t for t < to and u = to for t > to 
where to is a maximum integration time. The random integrand 
in the inequality is related to the random signal excess at the 
time s for an integration time u. The relation is: 




(10) (1/u) [R(s)/Rm(u) ] * ds > 1 




(11) 10 log [R(s)/Rm(u) ] * = 2[SE(s;u) + X(s) ] 
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where SE(s;u) is the expected value of the signal excess at a 
time s for an integration time u and X(s) is the random 

component of the signal excess at the time s. In the model, 

X(s) is determined by a lambda-sigma jump process and SE(s;u) 
is determined by an expected value sonar equation with a 
recognition differential RD(u) = 10 log rju(u) . Here, rm(t) is 
the value of the signal-to-noise ratio that gives a probability 
of detection equal to . 5 for an integration time t and a 
specified probability of false alarm pf. With the signal 
detection process described by a Case II signal detection model, 
the detection index necessary to give the required operating 
point (Pf/.5) is related to the integration time t and the 
signal-to-noise ratio rnv(t) by: 

(12) d = u-(BW) [r„,(t)]* 

where BW is the bandwidth of the receiver. For a spectrum 
analyzer, BW would be the bandwidth corresponding to a given 
frequency resolution and d would be the detection index 
required in order to be at the operating point (pf,.5) for a 
signal that was contained within a bandwidth BW. Since d in 
Equation 12 must be the same for t = u and t = tQ, 

(13) RD(u) = 5 log(to/u) + RD(to) 

where tg is the maximum integration time. Then, since 
SE(s;u) - SE(s;to) = RD(tg) - RD(u) , by using Equation 13 and 
Equation 11, Relation 10 becomes: 

ft (1/5) [X(s) + SE(s;to) - 5 log(to)] 

(14) Jt-u 10 ds > 1 
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where as above the time origin is chosen so that t > 0, the 
integration time u = t for t < to and u = to for t > to 
and where SE(s;to) is the expected value of the signal excess 
at the time s for a recognition differential determined by an 
integration time to* In an encounter, detection occurs the 
first time that Relation 14 is satisfied. 

As is pointed out in Reference 11, the appeal of the Third 
Passive Sonar Encounter Detection Model relative to the Second 
and First Passive Sonar Encounter Detection Models is that it 
appears to more closely describe the detection process in passive 
sonar detection systems that display their processed data to an 
operator in a continuous manner over a time window of duration 
tQ. However, results reported in Reference 12 indicate that the 
difference between the three models may not be significant in 
some types of encounters. 
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VIII. straight Line Encounters 

Suppose a target's detectability depends on its range 
from a detection system and that the probability of detection is 
effectively zero beyond a range r^^ for any target azimuth. In 
this report, an encounter between the target and the detection 
system is the event that the range between the target and the 
detection system is less than or equal to In addition, 

suppose r^ is small enough so that when the target and the 
detection system are having an encounter they can be considered 
to be moving on planes parallel to a tangent plane to the earth's 
surface at some point in their vicinity. If this is the case, 
then while the target and detection system maintain a constant 
course and speed during an encounter, the encounter is called a 
straight line encounter. 

A straight line encounter can be described in terms of a two 
dimensional rectangular coordinate system whose plane is parallel 
to the tangent plane to the earth. If the coordinate system is 
stationary relative to the detection system with the detection 
system located at the origin and is oriented so that the target's 
motion is parallel to the y-axis and is in the positive 
y-direction, then the target's x-coordinate during a straight 
line encounter will be constant. The constant is equal to the 
target's horizontal range at the closest point of approach (CPA) 
on the straight line track on which the target is moving relative 
to the detection system during the encounter. This range is 
called the target's lateral range. 



30 



A complete straight line encounter is a straight line 
encounter that begins at a range from a detection system that is 
greater than or equal to rj„ and continues past CPA to a range 
from the detection system that is again equal to or greater than 
rj„. Let p(x) be the cumulative probability that a target is 
detected by a detection system in a complete straight line 
encounter in which the target's lateral range is x. Then the 
function p(x) defines what is called a lateral range curve or 
lateral range function. 

Let p be the probability that a target is detected during a 
complete straight line encounter. If the lateral range of a 
target in a straight line encounter is assumed to be a continuous 
random variable X with a uniform distribution with fx(^) “ 
for |x| < a/2 and p(x) = 0 for |x| > a/2, then the 
probability that a target will be detected during a complete 
straight line encounter is given by; 



where the limits of integration can be used since the value of 
p(x) is zero for |x| > a/2. Equation 15 suggests a measure of 
a detection system's capability to detect a target in a straight 
line encounter. The measure W is called sweep width and 



(15) 




dx 



(16) 




dx. 



With this definition. Equation 16 becomes: p = (1/a) »W. 
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IX. Two Intermittent Signal Encounter Models 

In the intermittent signal encounter models that are 
described in this section, an encounter is a complete straight 
line encounter, and during an encounter a target either emits a 
signal (an acoustic transient) or its presence (a visible 
submarine mast) is the cause of a signal at various times. Two 
cases are considered: In the first case, the signals occur 

periodically, the signals are of length St and the time between 
the occurrence of signals is t where t > St. In the second 
case, St = 0 (the signals are instantaneous) and the signals 
occur at times determined by a Poisson process for which the 
expected time between signals is equal to T. In the model, the 
detectability of a target signal depends on a target's horizontal 
range from a detection system, but on no other factors. If a 
signal occurs while the target is within a range r, it will be 
detected. For a continuous signal, the lateral range function of 
a detection system for a target is: p(x) = 1 for |x| < r and 

p(x) =0 for |x| > r where the horizontal range r is 
determined by the characteristics of the detection system and the 
target. The geometry for an encounter is shown in Figure 3 
below. 

For intermittent signals, the length of a target's track 
relative to a detection system on which a signal will be detected 
is 2*(r* - x*)5 + w*5t where w is the speed of the target 
relative to the detection system. So, a target's exposure time 
during an encounter is (2/w) • (r* - x*)^ + St. 
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For periodic signals, there are two cases. In the first 
case, r > w(t - 6t)/2. In this case, the signals result in the 
following lateral range function: 
p(x) = 0 for |x| > r 

(17) p(x) = 1 for |x| < {r* - [w(t - St)/2]*)^ 
p(x) = [2/ (w* t) ] • ( r* - x*)^ + 6t/T otherwise 

In the second case, r < w*(t - 5t)/2 and the middle equality 
in Equation 17 does not apply. 

For signals that are instantaneous and whose occurrence is 
determined by a Poisson process, the signals result in the 
following lateral range function: 

p(x) = 1 - exp{ -[2/ (w T) ] • (r* - x*)») for |x| < r 

(18) 

p(x) = 0 for |x| > r. 




Figure 3 . The encounter geometry for the two intermittent signal 
models described here. 
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For signals whose occurrence is determined by a Poisson 
process and for which St > 0, signals can overlap. If this is 
allowed, then Equation 18 can be modified to describe this case 
by adding St/i to the term in the exponent of Equation 18 that 
is within the square brackets. In particular, note that this 
modified Equation 18 can be approximated by the bottom equality 
in Equation 17 when (2/WT)(r* - x*)^ + St/T « 1. This 
implies that when the expected time t between signals is large 
relative to the exposure time (2/w) (r* - x*)^ + St, the periodic 
signal model and the Poisson random signal model are effectively 
equivalent. 

If St = 0 and r < (wt)/2, for the periodic intermittent 
signal, W = Tr*r*/(VT). With the above approximation, this is 
also the sweep width for the Poisson random signal model. 
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X. A Random Search Model 



A search of a region in which a target moves on a track 
that consists of a number of straight segments placed in such a 
way that in a limiting sense every section of the region is 
equally likely to be searched on a segment is referred to as a 
random search in Reference 6. Representation of a search region 
with the track segments that could be imagined to be the tracks 
of a random search are shown in Figure 4 . 




Figure 4 . A search region and a track that could be described as 
a random search track. 

Two developments of a model to describe this kind of search 
are contained in this section. The first development is based on 
the following conditions: 1. A target is at a fixed position 

within a defined search region. 2. A searcher's track is a 
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sequence of straight line segments that are within the search 
region. 3. The searcher's detection system is such that while on 
a track segment, a rectangle is searched that is contained within 
the search region, is of length equal to the length of the track 
segment and is oriented so that its long axis is parallel to the 
track segment. 4. The probability that the searcher's detection 
system will detect a target while on a track segment with a 
search rectangle that does not contain the target is zero. The 
probability that the searcher's detection system will detect a 
target while on a track segment with a search rectangle that 
contains the target is p(x) where x is the target's lateral 
range for the track segment and p(x) is the lateral range curve 
for a complete straight line encounter lateral range x. A 
representation of a search rectangle is shown in Figure 5 below. 
5. The track segments are located in such a way that the event 
that the target is within the search rectangle associated with a 
track segment is independent of the event that the target in the 
search rectangle associated with any other track segment. And 
the probability of the event is equal to the ratio of the area of 
the search rectangle to the area of the search region and, given 
a target is within a search rectangle, its position is uniformly 
distributed over the rectangle. 

Condition 4 implies that the random search model is based 
on the concept of a complete straight line encounter. The 
definition of an encounter that is intended here is that given in 
Section VI. This implies that in the random search model the 
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time to resolve a false alarm is zero. However, for the model, 

P(j and pf are considered to be determined by some criterion 
such that pf is less than one. Consequently, although the time 
to resolve a false alarm is zero in the model, the cost 
associated with a false alarm is not zero. (A simple model that 
accounts for the time to resolve false alarms is described in 
Reference 13.) Condition 4 also implies that when a searcher is 
on a track segment with a search rectangle that contains a 
target, the encounter is a complete straight line encounter. And 
Condition 5, which can be considered to specify a random 
arrangement of the track segments, implies that when this is the 
case, for the complete straight line encounter, the target's 
lateral range is a random variable that is uniformly distributed 
between -b/2 and b/2 where b is width of the search 
rectangle (the dimension of the rectangle perpendicular to the 
associated track segment) . 



Swept Area 



Track 

Figure 5. A track segment and its associated search rectangle, 
that could correspond to a search with an aircraft mounted 
infrared detection system. 
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Based on the above considerations, the probability that a 
target will be detected while a searcher is on a track segment 
with an associated search rectangle that contains the target is 
given by: 



where fx(^) ~ for -b/2 < x < b/2 and fx(^) ~ ® 

p(x) = 0 otherwise. Note that the left side of Equation 19 
applies to any complete straight line encounter in which the 
target's lateral range for the encounter is considered to be a 
random variable with a distribution determined by the probability 
density function . If it is not given that the target is 

within the search rectangle associated with a track segment, then 
the unconditional probability that the target will be detected on 
the track segment is given by: (W/b)*(5A/A) where 5A is the 

area of the search rectangle associated with the track segment 
and A is the area of the search region. With 1 the length of 
the rectangle, 5A = b*l and the probability becomes: (W*l)/A. 

Then, since the event that the target will be in the search 
rectangle of a track segment is independent of the event that it 
will be in the search rectangle of any other track segment, the 
probability p that a random search consisting of m track 
segments will detect the target is given by: 

1 - [1 - (W*l 3 ^)/A][l - (W*l 2 )/A] ••• [1 - (W*lj^)/A] where 1^ 
is the length of the i^^ track segment. The probability is 
also given by: p = 1 - exp{2 ln[l - (W«li)/A]) where the sum 



(19) 




fx(x) dx = W/b 
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• • • 



index i = 1 , 2 , ••• ,n. If (W‘li)/A « i for i = 1 , 2 , ••• 

then this expression can be approximated by: 

(20) p = 1 - exp[-(W‘l)/A] 

where 1 = S is the track length of the search. Equation 20 

is known as the random search formula. 

The second development of the random search formula is based 
on Equation 5 and a detection rate for a random search given by: 

T (t) = W‘V(t)/A. With this detection rate and Equation 5 , the 
random search formula is given by: 

(20a) P(T < t) = 1 - exp {-(W‘l (t) ]/A) 

where l(t) is the track length for a random search that starts 
at time 0 and ends at time t and 



(20b) l(t) = Jo v(s) ds. 

Replacing P(T < t) by p and l(t) by 1 gives Equation 20. 
In the form of Equation 20a, the random search formula indicates 
explicitly the relation between the probability of detection and 
the duration of a random search. Note that Equation 20a implies 
that the sweep width is independent of speed over the range of 
speeds in the encounter. 

As an example application, consider the periodic signal 
model of Section IX with 5t = 0, r < (wt)/ 2 and v(t) = v. 
For this case, p = 1 - exp[-(Ti‘rVA) • (t/t) ] . 

Reference 14 contains an example of an application of the 
method used in the second development of the random search 
formula to a random search where the search region expands with 
time. 
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XI. Ladder and Barrier Search Models 

In some barrier searches, the barrier search track is a 
ladder search track relative to a reference system that moves 
with the target. This fact is used in the barrier search model 
development that follows the two ladder search model developments 
below. The first ladder search model is referred to as an ideal 
ladder search model because of the idealizations that are 
involved in its description of a ladder search. The second 
ladder search model is referred to as a degraded ladder search. 

It can be considered to describe a ladder search track in which 
navigational errors result in omissions and overlaps in coverage. 

An Ideal Ladder Search Model: The model is based on the 

following conditions; 1. A ladder search region is a rectangle 
that contains a fixed target. 2. During a search of the region, 
the searcher's detection system moves on a set of m parallel 
track segments of length b separated by a distance s. 3. As 
the detection system moves along a track segment, it searches a 
rectangular strip of length b and width s within the search 
region. 4. The m rectangular strips that correspond to the m 
track segments completely cover the ladder search region with no 
overlap. 5. If a target is within the rectangular strip 
corresponding to a track segment, then there will be a complete 
straight line encounter between the target and the detection 
system when the detection system moves along the track segment 
and the lateral range of the encounter will be uniformly 
distributed across the width of the strip. If the target is not 
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in the rectangular strip, then there will not be an encounter and 
the probability that the target will be detected while the 
detection system is on the track segment is zero. 




Figure 6. A schematic representation of a ladder search geometry 
for a case in which the ladder search track segments are 
superimposed on and bisect their corresponding rectangular 
strips . 

Since targets outside of the rectangular strip that 
corresponds to a track segment cannot be detected while a 
detection system is on the track segment because of Condition 5, 
in the model, the sweep width W of a searcher's detection 
system must satisfy the relation W < s . In particular, W = s 
only holds when the detection system detects a target that is in 
a rectangular strip with probability one for any target lateral 
range. This kind of detection system is sometimes referred to as 
a cookie cutter detection system. However, this terminology can 
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be misleading since it suggests the detection system detects 
equally well for all azimuths. But this is not a requirement on 
the system in order that W = s. 

The ideal ladder search model implies that if the conditions 
of the model are satisfied, then the probability p that a 
target will be detected by a an ideal ladder search is given by: 

(21) p = W/s 

where W/s < 1. The quantity W/s is called the coverage 
factor in this case. 

A Degraded Ladder Search Model: The above model implies 

perfect navigation in addition to other idealizations. A model 
of a ladder search is given in Reference 6 that could be used 
for cases in which this is a poor assumption. The model which is 
referred to here as a degraded ladder search model can be 
considered to describe navigational inaccuracies in terms of 
omissions and overlaps of the rectangular strips. It can be 
developed as follows: Consider a random search in the ladder 

search region whose track length is equal to the search track 
length required to complete an ideal ladder search, that is, a 
track length 1 = m*b. The degraded ladder search model 
describes the result of omissions and overlaps in a ladder search 
to be such that the probability of detection for this random 
search is equal to the probability of detection for the degraded 
ladder search. Consequently, since the area of the ladder search 
region is m*s*b, for the degraded ladder search model: 

(22) p = 1 - exp(-W/s). 
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Here, the requirement that the coverage factor W/s < 1 for 
Equation 21 can be relaxed. However, it should still be 
considered as an approximate condition. 

The condition that the target be fixed within the 
rectangular search region is critical to both Equation 21 and 
Equation 22. However, these results are also applicable to a 
search for a moving target under the conditions that are 
described next. 

A Barrier Search Model: A target moves with a constant 

course and a constant speed u. Both the target's course and the 
target's speed are known by a searcher. The searcher establishes 
a barrier of width b that is perpendicular to the target's 
track and moves on the barrier with a speed v > u. The barrier 
is designed so that in a reference system relative to the target 
the barrier search is a ladder search that satisfies the 
conditions for a ladder search that are given above. There are 
two cases to consider; 1. The barrier is established in front of 
the target. 2. The barrier is established behind the target. 

From the search geometry for a barrier established in front 
of the target, it can be seen from Figure 7 below that 
8 = sin“^(u/v) and d = vt where t = s/(v + u) is the time to 

move from one search leg to the next. The angle 0 and the 

perpendicular distance d that depend on u, v and s, and the 

width of the barrier b are the quantities that are required in 
order to establish the barrier operationally. 
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Target 




b 

Figure 7. A barrier search track shown for a barrier established 
in front of the target. The track is shown in a reference system 
fixed relative to the earth. 

For a barrier that is established in front of a target, one 
of three barrier types will result. A barrier's type is 
determined by the relation of the distance d to the distance 

g = ut where the time t = b/(v^ - u^ ) ^ is the time to complete 

a search leg (cross the barrier) . The barrier type is determined 
as follows: 1. For g < d, the barrier is an advancing barrier. 

2. For g = d, the barrier is a stationary barrier. 

3. For g > d, the barrier is a retreating barrier. 

For a barrier established behind the target, there is only 
one barrier type and it is called an overtaking barrier. For an 
overtaking barrier, 9 = sin~^(u/v) as for a barrier established 
in front of the target. But, for an overtaking barrier, 

T = s/(v - u) and d = vs/(v - u) . 
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Given the target crosses the barrier, the probability of 
detection for an ideal barrier search is given by Equation 21 
and the probability for a degraded barrier search is given by 
Equation 22 where the terminology refers to the nature of the 
ladder search in the reference system moving with the target. 

A discussion of an application of these two equations to a search 
for a magnetic anomaly target is given in Reference 15. 
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XII. A Target State Estimation Procedure 

A target state estimation procedure based on bearing 
observations is developed in this section that generates point 
estimates of a target's position and velocity vector coordinates 
in a rectangular coordinate system. The procedure is based on a 
model in which bearing errors are unknown and are not determined 
by random variables with known distributions. Because of this, 
confidence regions for the estimates are not generated by the 
procedure. However, for a moving target, it illustrates general 
characteristics of bearings only target motion analysis (TMA) . 

The model is defined as follows: 1. The target moves in a plane 

with a constant but unknown course and speed. 2. Observations of 
the target are made from known positions at known times. 3. The 
observations provide only target bearings with unknown errors. 

The model geometry is shown in Figure 8 . 



North 




observed G-i 
bearing 



[Xo(i) ,yo(i) ] observer 



X East 



Figure 8. The geometry of the target motion analysis model. 
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The procedure criterion is: For observations from n positions, 

choose target position estimates and target velocity component 
estimates and Uy that make the sum of the squares of the 

algebraic distance between the estimated positions and their 
corresponding observed bearing lines a minimum. From Figure 8, 
it can be seen that the algebraic distance can be written as 
‘^i = [^t(i) - Xo(i)]*cos - [Yt(i) - yo(i)]*sin 0^. Because of 
the requirement that the target move with constant course and 
speed during the encounter, the number of independent estimates 
is reduced from 2n to 4, , Uy and any two position 

estimates x^Cj), yt(j)* in the following development, j = 1 
and with i = 2,3, ,n the remaining estimates are given by: 

xt(i) = xt(l) + Ux*(ti - ti) and yt(i) = yt(l) + '^y* (ti - t^) . 
To determine "best" estimates of the target state parameters, 
take the partial derivative of the sum S = Z (d^) ' with respect 
to each of them. Then set the four partial derivatives equal to 
zero. This creates four linear equations in x^Cl)/ yt(l)/ 
and Uy whose solution are the desired estimates x^(l), yt(l)/ 
Ujj and Uy. In matrix notation, the equations can be 
represented by AX = B where the elements of X are: 
xii = xt(l), X 21 = yt(i)f ^31 = ^41 = ^ necessary 

condition for a unique solution for X is that n > 4. 

Otherwise, the determinant of A will be equal to zero. The 
procedure can also be used if a target's course and speed are 
constant and known and, in particular, if the target is 
stationary so that Uj^ and Uy both equal zero. In this case. 
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since the number of unknowns is two, the number of linear 
equations is also two and a necessary condition for a unique 
solution is n > 2. 

Now, suppose the observations are at positions and times 
that correspond to the positions and times of an observer moving 
on some constant course at some constant speed (including zero 
speed) . In this case, the observation position coordinates are 
related by the following equations; Xq( 1) = Xo(l) + Vjj(ti - t^^) 
and yo(i) = Yo(^) ^y('*^i “ '*^l) where Vj^ and Vy are the 

required velocity components of the observer. Using these 
equations of motion, the matrix equation AX = B can be 
transformed to the matrix equation Ax'= 0 where the elements 
of the matrix X are related to the elements of the matrix x' 
by the equations; = Xt(l) - Xo(l), X 21 = yt(l) ~ yo(^)/ 

^31 = ^x " ^x X 41 = ^y ■ 

Since the linear equations represented by AX' = 0 are 
homogenous, they do not have unique solutions and consequently 
neither do the equations represented by AX = B. However, if 
there is at least one observation whose time and position is not 
determined by the above equations of motion, then the 
transformation from X to X* cannot be made, and in general a 
unique solution for X can be found. If the observations are 
made from a platform that is moving with a constant course and 
speed, this condition can be achieved by either changing the 
course, the speed or both prior to completing the observations. 

Estimation models that describe bearing error as a random 
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variable provide a basis for deteinnining confidence regions for 
point estimates. A model is developed in Reference 16 that does 
this for either target bearing observations made from two or more 
points simultaneously or for a target that is stationary relative 
to the observation points. 
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XIII. Position Distributions That Change with Motion 



Target motion models provide a basis for determining 



position distributions that change with target motion. In this 
section, two classes of target motion models are considered. In 
the first class, a target moves in a plane with a constant course 
and speed and the course and speed are independent of the 
target's position. In the second class, a target moves in a 
plane but its course or speed changes during the motion. Three 
members of the first class are developed first. This is followed 
by a brief discussion of some models of the second class. 

Motion Models of the First Class: For the first class of 

motion models, the joint density function of the distribution 
that determines a target's coordinates X(t) and Y(t) at some 
time t > 0 can be determined by: 



where V = U^, W = Uy and Ujj and Uy are the random 
variables that determine the target's velocity components Uj^ 
and Uy and v = u^, w = Uy, q = x- vt and s = y - wt. 
Equation 23 can be developed as follows: To first order, 

fx(t) ,Y(t) Sx5y is the probability that a target's 

coordinates are in an element of area 6x6y and for given values 
of V, w and t, target positions in an element of area 6q6s 
will be translated to an element of area SxSy that is identical 
in form and size to 5q*5s. And so, to first order, 

^X(O) /Y(0) Sx5y fy^Yj(u,v) Sv5w is the probability that 

the target's coordinates at time 0 are in an element of area 
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Sq&s that is located such that the target's coordinates will be 
in the element of area SxSy at time t since x = q + vt and 
y = s + w*t. And, to first order, the sum of such probabilities 
for all pairs of values of v and w is also the probability 
that the target's coordinates at time t are in the element of 
area 6x6y. In the limit after equating the two expressions for 
this probability and cancelling the common factor 6x»5y, 

Equation 23 results. 

The First Motion Model: In the first model, X(0) and 

Y(0) are both independent normal random variables with means 
and jiy and equal standard deviations a. However, and Uy 

are not normal and they are not independent random variables. In 
this model, Ujj = u*sin # and Uy = u*cos i where § is the 
random variable that determines the target's course and u is 
the targets speed which is known. So, only a value for the 
random variable $ is required to determine the target's 
velocity. In the model, $ has a uniform distribution over the 
interval 0 to 2tt and it is convenient to chose the 
rectangular coordinate system so that the means H-y 

each equal to 0. Then, with the circular normal distribution 
determining the random position coordinates and with the 
distribution that is described above determining the random 
velocity components , in the coordinates u and 0, the integral 
of Equation 23 is a single integral over <p and the integrand 
of the integral is (l/2Tro*) exp[-(q* + s*)/2a*] (l/2ir) where 
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now q = X - u-t*sin <f> and s = y - u*t*cos <p. Integration 
gives fx (t) , Y (t) V ' t) as^ 

(24) (l/2Tia*) exp{-[x* + y* + (u*t)*]/2a*) IqCCx* + y*)^*u*t/a*] 

where t > 0 and Iq indicates the hyperbolic Bessel function of 
zeroth order. In Reference 6, fx(t) , Y(t) plotted for 

several values of t in terms of r=(x*+y*)^, the target's 
range from the origin. The plots show a characteristic of the 
distribution that can be indicated as follows; First, replace 
(x* + y*)^ by r in fx(t) , Y (t) • Next, multiply and 
then divide fx(t) , Y (t) exp(-r‘U*t/a*) • This gives: 

(25) l/(2iTa*)exp{-[l/(2a*) ] (r - u* t) * }Io (r *u* t/o* ) exp (-r ‘U* t/a* ) 
where t > 0. As noted in Reference 17, Iq ( z) *exp(-z) is a 
slowly decreasing function that asymptotically approaches 
1/(2 'JTz)* as z increases. The consequence of this is that a 
plot of fx(t) , Y(t) against r for values of t greater 
than 4*a/u has the appearance of a normal density function. 

A target's random rectangular coordinates X(t) and Y(t) 
and its random bearing 6(t) and range R(t) from the origin 
are related by; X(t) = R(t) ‘sin 0(t) and Y(t) = R(t) *cos 0(t). 
Using these relations, ^x(t) ,Y(t) transformed to 

the joint density function fR(t) ,0(t) random 
variables R(t) and 0(t) . To do this, replace x* + y* by r* 
in Expression 24. Then multiply the resulting expression by r, 
the Jacobian of the transformation. This gives the expression: 

(26) (l/27r) (r/2Tra*) exp(-[r* + (u*t)*]/2a*) Io(r*u*t/a*) 
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where 0 < r and 0 < a < 2 tt. And, this is the joint density 
function fR(t) ,9(t) inspection, the marginal density 
function of ®(t) is 1/2'jr over the interval 0 

to 2ir and the marginal density function fR(t) R(t) 

is Expression 26 multiplied by 2ir. 

The Second Motion Model: In the second model, X(0) and 

Y(0) are independent normal random variables with means jijj and 
M.y and standard deviations Ojj and Oy that determine a 
target's random position coordinates at time 0. And Ujj and 
Uy are independent normal random variables with means Uj^ and 
Uy and standard deviation that determine a target's random 

velocity components. Because of these conditions, the target 
coordinates are X(t) = X(0) + Ujf*t and Y(t) = Y(0) + Uy*t at 
time t. This implies that X(t) and Y(t) are independent 
normal random variables with means + Ux*'^ M^y ^y*^ 

and with standard deviations Oy + and Oy + o^*t*. The 

model describes a bivariate normal position distribution whose 
center moves with a constant velocity determined by Uy and Uy 
and which becomes more and more circular as its standard 
deviations increase with the passage of time. Although the 
targets joint density can be found by using Equation 23, this 
procedure is more direct. For another discussion of the first 
and second models, see Reference 7. 

The Third Motion Model: In the third model, the target is 

at the origin of a rectangular coordinate system at time zero. 
After that, its position is uniformly distributed on a circular 
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disk of radius Uj„'t centered at the origin. This implies that 
the joint density function of the position distribution is: 

(27) fx(t) ,Y(t) (X/y't) = 

for t > 0 where x* + y* < u^*t’ and that the joint density 
function of the distribution of the random variables Q(t) and 
R(t) that determine a target's bearing from the origin is: 

(28) fR(t) , 0 (t) (J^/Ci't) = r/(TT*u^*t*) 

for t > 0 where 0 < r < Uj„*t and 0 < a < 2 * 7 T. 

Since the values of r and a are independent and the joint 
density function is equal to the product of l/( 2 *ir) and 
2*r/(u^*t*), the random variables R(t) and 0(t) are 
independent and fR(t) ('^) “ 2 »r/(u^*t*) where 0 < r < Ujf,*t and 
f 0 (t) (ci) = 1/(2*7 t) where 0 < a < 2'7T. These two marginal 
distributions can be achieved by choosing at time 0 a course <p 
from the uniform distribution with density: f$( 0 ) = l/( 2 *ir) 

where 0 < 0 < 2‘tt and a speed U from the triangular 
distribution with density: fu('^) = 2 *u/u^ where 0 < u < Uj„. 

These choices define the third motion model. 

Motion Models of the Second Class: For the second class of 

motion models, a target's course or speed or both can change. In 
general, a monte carlo simulation method is required in order to 
determine a position distribution that is based on such a model. 
As an example of cases in which the distribution can be described 
analytically, see Reference 18 and Reference 19. As an example 
of a case in which it can not, suppose a target's initial 
position is described in terms of a number assigned to a 
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subregion in the xy-plane where the number assigned represents 
the probability that the subregion contains the target at an 
initial time. In addition, suppose for each subregion a course 
and speed distribution is determined by assigning numbers to 
course and speed pairs where a number represents the probability 
the target will have the course and speed at the initial time 
given it is in the subregion at that time. Next suppose for each 
course and speed pair there is a time distribution that 
determines the duration of the course and speed pair and that the 
time distribution is determined by a number assigned to each 
discrete time point where the number represents the probability 
that the target's course and speed pair will be determined by a 
new course and speed distribution. By extending this kind of 
procedure and then implementing it in a monte carlo simulation, 
one can generate complex position distributions that describe a 
target's position at discrete time points. 
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XIV. Position Distributions That Change with Search 

For the models that are considered here, a target is 
within a region that has been divided into subregions or cells. 
And, for each cell, a number has been assigned to the cell that 
is interpreted as the probability that the target is within the 
cell. The set of these probabilities constitute a position 
distribution. Suppose information becomes available that a 
search has been conducted for the target and that the target has 
not been detected. Or suppose the information is that the target 
has been detected. In the first case, negative information is 
available that can be used to modify the position distribution. 

In the second case, positive information is available that can be 
used to modify the position distribution. 

Position Distributions and Negative Information: For a 

region that contains a target and consists of n subregions, let 
the event = {the target is in the i^^ subregion}. And let 

the event C = (no contact}. Then, given no contact in a search 
of the region, the targets 's position distribution can be 
modified as follows: 

(29) P(Si|c) = P(C|Si) •P(Si)/P(C) 

where i = 1,2, ••• ,n and P(C) = E P(clSj} •P(Sj) with the sum 
index j = 1,2, ••• ,n. Note that Equation 29 can be obtained 
by using Bayes theorem. To illustrate how Equation 26 might be 
used, suppose that a search in a subregion is considered to be a 
random search and that the sweep width of detection system 
against the target depends on subregion being searched. For this 
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case, let Aj^ be the area of the subregion and let Wj^ be 

the sweep width in that subregion. Then, given no contact in a 
search of a subregion, P(C|Sj^) = exp[-(W£*lj^)/Aj^] where Ij^ is 
the track length of the searcher in the i^^ subregion. Given 
values for P(Sj^) , W^/Aj^ and Ij^ for i = 1,2, ••• ,n, a 
position distribution can be determined that has been modified by 
the negative information. 

Position Distributions and Positive Information: In the 

case of positive information, the event C = {a contact) occurs. 
Then, given a contact in a search of the region, the target's 
position distribution can be modified as follows: 

(30) P(Si|C) = P(c|Si) -P(Si)/P(C) 

where P(C) = E P(c|Sj)*P(Sj) and the sum index j = 1,2, ••• ,n. 

To illustrate a way Equation 30 might be used, suppose that a 
contact is a line of bearing detection or an omnidirectional 
sensor detection and that the cells are range cells. Then, with 
rj^ a range identifying the i^^ cell and R a random variable, 
P(Sj^|C) becomes P(R = rj^lc), P(c|Sj^) becomes P(c1r = rj^) and 
P(Sj^) becomes P(R = rj^) . In a continuous analog of this, 
fj^(r|C) replaces P(Sj^|c), P(c|r = r) replaces P(c|Sj^) and 
fp(r) replaces P(Si) . 

Position Distributions and Uncertain Information: The event 

C is the union of the two events: T^ = (a true contact) and 

Fj, = (a false contact). And the event C is the union of two 
events: T^ = (a true no contact) and Fq = (a false no contact) 

The correspondence between C, C, Sj^, Sj^, T^, F<^, T^ and 
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Fq and the events of the Venn diagram in Figure 1 are as 
follows: C corresponds to D^, C corresponds to Dq, Sj^ 

corresponds to Sj^ corresponds to Hq, corresponds to 

n , Fq corresponds to (D]^ o Hq) , corresponds to 

(Dq n Hq) and F 3 corresponds to (Dq n H^) . 

Using the above relations, after a search of a region that 
has resulted in a contact, the target's position distribution can 
be defined by: 

(31) P(Si|C) = P(Si|Tc) -P(Tc|C) + P(SilFc) -PCFcIC) 

where i = 1,2, ••• ,n since = T,- n C, Fq = Fq n C and 

P(Si|c) = [P(Si n Tc) + P(Si n Fc) ]/P(C) . The probability 
p = P(Tq|C) has been called the credibility of the contact. In 
terms of p Equation 31 becomes: 

(32) P(Si|c) = P(Si|Tc)-p + P(SilFc)-(l - p). 

Again using the above relations and p = P(T 3 |c), after a 
search of a region that has resulted in no contact, the target's 
position distribution can be defined by; 

(33) P(Si|c) = P(Si|Tc)*p + P(Si|Fc)*(l - P) 
where i=l, 2 , ••• ,n. 

The values of P(Tc|Sj) and of P(FclSj) = 1 - P(Tc|Sj) 
are determined by the characteristics of the search in the 
subregion for j = 1,2, ••• ,n. With values for these two 

probabilities, P(Sj^|Tc) = P (T^ | Sj^) • P (Sj^) / [E P(Tc | S j ) • P(Sj ) ] and 
P(Si|Fc) = P(Fc|Si) -P(Si)/[S P(Fc|Sj) -P(Sj) ]. 

In one positive information model, p is determined 
subjectively based on the detection system and the nature of the 



58 



search and P(Si|Fc) = P(Sj^) . This choice for P(Si|Fc) could 
be based on the argument: Given the detection system and the 

nature of the search, for a search that ends in a false contact, 
the target was not yet detectable. Consequently, a search that 
ends in a false contact supplies no information about a target's 
location. This argument implies that there have been no missed 
detections during the search. 

In a negative information model, in keeping with the above 
choice for P(Sj^|Fc), one might choose P(Sj^|Tc) = P(S^) . This 
choice for P(Sj^|Tc) could be based on an argument that is a 
parallel to the one for the choice above for PCS^If^): Given 

the detection system and the nature of the search, for a search 
that ends in a true no contact, the target was not yet 
detectable. Consequently, a search that ends in a true no 
contact supplies no information about the target's location. In 
keeping with the above comment, this argument implies that if 
false alarms occurred during the search, then they were resolved. 

If p = 1, then C = T^ (a detection) and, using the above 
expression for P(Sj^|T£,), Equation 32 becomes identical to 
Equation 30. Note that 1 - p corresponds to P(Ho|Di) and not 
to P(Dj^|Hq) which is pf. 

If p = 0, then C = Fc (a missed detection) and, using the 
above expression for PCSj^lF^), Equation 33 becomes identical 
to Equation 29. 
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XV. Search Models and Search Theory 

Search theory provides a basis for determining optimal 
search plans for a target whose motion and location are 
determined within some bounds. Here, an optimal search plan is 
one for which the probability of finding a target within a given 
length of time is a maximum, the expected time to find a target 
is a minimum given the target is found or a search plan for which 
some other optimal search criterion is satisfied. 

Search theory results are based on models of the search 
process. To the degree that a search model describes a search 
process, an optimal search plan for a target that is based on the 
search model should provide guidance for the development of an 
operationally feasible search plan. However, because of the 
limitations of analytical search models, an optimal search plan 
that is based on an analytical search model may give only initial 
guidance in this regard. The optimal search plans that are 
described below illustrate this. The search plans are based on 
the random search model. Because of this, the requirement on the 
location of search track segments is not realizable and the time 
to resolve false alarms is ignored. 

Optimal search plans based on search models implemented 
through a monte carlo simulation are not considered here. 

However, with sufficient information, such plans have the 
potential of being both implementable and more optimal in a real 
sense than an optimal search plan based on an analytical search 
model . 
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Three Ophinal Search Plans: The three optimal search plans 

differ in their criterion for an optimal search plan. However, 
each one is based on the following search model: A target is 

fixed at some point in a region that consists of n subregions. 

A search in a subregion is a random search in the sense of the 

definition in Section X and a searchers sweep width there is a 

constant. In addition, a search of a subregion will not detect a 
target which is in another subregion. To determine a plan, let 
Si = {the target is in subregion i} for i = 1,2, ••• ,n and 
let Pi = P(Si) be the prior probability that the target is in 
the i^*^ subregion. Let Wi be the sweep width in the i^*' 

subregion. Let &i = Ai/Wi where Ai is the area of the i^^ 

subregion and 6i is the expected track length to find the 
target by a search of the i^^ subregion given the target is in 
the i^^ subregion, a characteristic length. The probability P 
that the target will be detected by a random search is given by: 

(34) P = S [1 - exp(-li/Si) ] -Pi 

where the sum index i = 1,2, ••• ,n and Ij^ is the track 
length of the search in the i^^ subregion. 

The first criterion: Choose 1^ so that P is a maximum 

subject to the two constraints: 1. 1 = S Ij^ and 2. 1^ > 0 

where the index i = 1,2, ,n. Determining this choice is a 

nonlinear optimization problem whose solution is given in 
Reference 20. It is: 

l^/6i = ln(pi/6i) - L(k) i = 1,2, ••• , k 

(35) 

l|/5i =0 1 = k+l,k+2, ,n 
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where L(k) = (1/Z6 j ) *2[ 6 j • In (pj/S j ) ] - 1/S6j and the sum index 
j = 1 , 2 , ••• ,k, where the subregions are relabeled so that the 
following order relation holds: Pi/S^ > P 2 /S 2 > ••• > Pn/^n 

where k is chosen so that for k+i the solution for Ik+i 
using L(k+1) is either negative or zero. 

The second criterion: Choose 1^ so that P is a maximum 

subject to the two constraints: 1 . c = 2 c^ and 2 . Cj[ > 0 

where the index i = 1 , 2 , • ♦ • , n, cj^ = kj^ • 1 ^ is the cost of the 
search in the i^^ subregion and k^ is the cost per unit track 
length in that subregion. For this criterion, the solution to 
the corresponding nonlinear optimization problem can be obtained 
from Equation 35 by replacing 5^ by = kj[*5i and labeling 
the subregions so that Pi/^i > P 2 /^i ^ *** ^ Pn/^i* basis 

for this can be seen by replacing li/Si by its equivalent 
Ci/ei in the exponential term in Equation 32. 

The third criterion: Choose 1 ^ so that the expected 

utility of the search is a maximum subject to the two 
constraints: 1 . 1 = 2 Ij^ and 2 . Ij^ > 0 where the index 
i=l,2, ••• ,n. For this criterion, the solution to the 
corresponding nonlinear optimization problem can be obtained from 
Equation 31 by replacing p^ by q^ where q^ = Uj^»p^ and u^ 
is the utility of detecting the target given it is in the i^^ 
subregion. And, in addition, labeling the subregions so that 
qi /61 > q 2/&2 > **• > qn/^n* basis for this can be seen by 

multiplying the summation term in Equation 34 by Uj[ so that 
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the resulting equation gives the expected utility of the search 
given the utility of not detecting the target is zero. 

Equation 35 can be used to determine an order of search for 
the subregions which will effectively minimize the expected track 
length required to detect a target given it is detected. To do 
this, divide the available track length 1 into units small 
enough so that with a single unit only the subregion would be 
searched. Then allocate one unit to the search of the 
subregion. If the search is unsuccessful, determine the optimum 
allocation for two units. Then search with a second unit so that 
the first search with the first unit plus the second search with 
the second unit satisfy the optimum allocation for two units. If 
the search is unsuccessful, continue in this fashion until either 
the target is found or all the track length is expended. That 
this allocation order will effectively minimize the expected 
track length required to detect a target given it is detected can 
be argued as follows; Let L be the track length at detection, 
let lu be a unit of track length and let n be the number of 
units. Then the value of the probability P(L < i*lu) that the 
target will be detected on or before the i^^ step of the search 
for the given allocation order will be greater than or equal to 
its value for any other allocation order with the same allocation 
step size. Since the value of P(L < 1) will be equal to its 
value for any other allocation order of the optimum allocation 
and since P(L < i*lulL < 1) = P(L < i*lu)/P(L < 1), the value of 
the distribution function FL(i*ly|L < 1) = P(L < i^l^lL < 1) 
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will be greater than or equal to its value for any other 
allocation order. This implies that the expected track length 
given detection E(L|l < 1) = 2 [1 - FL(i«lu|L < 1)] where the 
sum index i = 1,2, ,n is effectively a minimum for the 

given allocation order. A search based on the optimum allocatic 
given by Equation 35 and the given allocation order is 
equivalent to the following search: After an allocation of trac 

length ly and an unsuccessful search, new values for P(Sj^) ai 
calculated using Equation 29 and then Equation 35 is used witt 
these new values to determine the next optimum allocation. A 
discussion of this procedure is given in Reference 6. And an 
example of its application is given in Reference 21. 

Equation 35 also defines an optimal search plan for a 
detection system that searches beams and can be described by 
Equation 33 by replacing 1^ by t^ where tj^ is the time tt 
i^h beam is searched and by replacing Sj^ by Tj^ where T|, 
characteristic time, is the expected time to detect the target I 
a search of the i^^ beam given the target is in the i^^ beam. 

For a more extensive discussion of search theory and its 
application to military operations research, see Reference 22. 
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